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Abstract 

In the framework of classical Glauber approach the analytical expressions for the 
variance of the number of wounded nucleons and binary collisions in AA interac- 
tions at given centrality are presented. Along with the optical approximation term 
they contain the additional contact terms, arising only in the case of nucleus-nucleus 
collisions. The magnitude of the additional contributions, e.g. for PbPb collisions 
at SPS energies, at some values of the impact parameter is larger than the con- 
tribution of the optical approximation, with their sum being in a good agreement 
with the results of independent Monte-Carlo simulations of this process. Due to 
these additional terms the variance of the total number of participants for periph- 
eral PbPb collisions and the variance of the number of collisions at all values of the 
impact parameter exceed several times the Poisson ones. The correlator between 
the numbers of participants in colliding nuclei at fixed centrality is also analytically 
calculated. 



1 Introduction 

At present the considerable attention is devoted to the experimental and theoretical inves- 
tigations of the multiplicity and transverse momentum fluctuations of charged particles 
in high energy AA collisions (see and references therein). One expects the increase 

of the fluctuations in the case of freeze-out close to the QCD critical endpoint of the 
quark-gluon plasma - hadronic matter phase boundary line [HI [9j. 

The aim of the present paper is to draw an attention to another factor leading to the 
increase of the fluctuations in the case of AA interactions. Namely the increase of the 
fluctuations of the number of participants and binary collisions due to multiple contact 
nucleon interactions in nucleus-nucleus collisions. 

Clear that these fluctuations lead to fluctuations in the number of particle sources and 
so directly impact on the multiplicity and transverse momentum fluctuations of produced 
charged particles and also on the correlations between them (see, for example, |10j-|17j). 

In the present paper the analytical expressions for the variance of the number of 
wounded nucleons and binary collisions in given centrality AA interactions are obtained 
taking into account the multiple contact NN interactions (so-called loop contributions). 
The calculations are fulfllled in the framework of classical Glauber approach [18], having 
a simple probabilistic interpretation [191 120] • In contrast with purely Monte-Carlo simu- 
lations the analytical calculations enable to understand the origin of increased values of 
the fluctuations. 



1 



As a result we demonstrate that the muhiple contact NN interactions in AA scattering 
lead in particular to the fact that, e.g. for PbPb collisions at SPS energies, the variance of 
the total number of participants for peripheral collisions and the variance of the number 
of collisions at all values of the impact parameter exceed a few times the Poisson ones. 

The paper is organized as follows. In section 2 in the framework of classical Glauber 
approach we present the analytical expression for the variance of the number of wounded 
nucleons in one of the colliding nucleus at a fixed value of the impact parameter. Along 
with the well known optical contribution (which depends only on the total inelastic NN 
cross-section) in the case of nucleus-nucleus collisions there is the additional contact term, 
depending on the profile of the NN interaction probability in the impact parameter plane. 

In section 3 we calculate the correlator between the numbers of participants in colliding 
nuclei at fixed centrality and as a consequence find the variance of the total (in both nuclei) 
number of participants. 

In section 4 in the framework of the same approach we present the analytical expression 
for the variance of the number of NN binary collisions in given centrality AA interactions. 
Along with the optical approximation term it also contains other terms, which occur the 
dominant ones. These terms also correspond to the multinucleon contact interactions and 
arise only in the case of nucleus-nucleus collisions. 

The derivations of all formulas are taken into the appendices A, B and C. 

All over the paper the results of numerical calculations are presented with the purpose 
to illustrate the obtained analytical results. We control also the results of our analytical 
calculations comparing them with the results obtained by purely Monte-Carlo simulations 
of the nucleus-nucleus scattering. 

Note that we restrict our consideration by the region of the impact parameter /5 < 
Ra + Rb, where the probability of inelastic interaction c^^W) nuclei with radii 

Ra and Rb is close to unity. 

2 Variance of the participants number in one nucleus 

At first we consider the variance \^[A^^(/3)] of the number of participants N^{f3) (wounded 
nucleons) at a fixed value of the impact parameter /3 in one of the colliding nuclei A. 
In the framework of pure classical, probabilistic approach to nucleus-nucleus collisions, 
formulated in [IB], we find for the mean value and for the variance of N^[I3) the following 
expressions (see appendix 1X1): 

{N^m = AP{(5) , (1) 

V[N^m = AP{P)QiP) + AiA - l)[Q(^2)(/3) - Q'm , (2) 

where P{(3) = 1 — Q{(3). For Q{(3) and Q^^'^\(3) we have (all integrations imply the 
integration over two-dimensional vectors in the impact parameter plane): 

Q{(3) = j da TA{a)[l - /B(a + /3)]^ , (3) 

Q(^2)(^) = j da,da2TA{a,)TA{a2)[l - /B(ai + /3) - fB{a2 + P)+9B{ai + P,a2 + P)f (4) 
with 

fsia) = J dbTB{b)a{a-b) , (5) 
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gB{ai,ci2)= / db TB{b)a{ai-b)a{a2-b) 



(6) 



Here Ta and Tb are the profile functions of tlie colliding nuclei A and B. The a{a) is 
the probability of inelastic interaction of two nucleons at the impact parameter a. We'll 
imply that o"(a), Ta and Tb depend only on the magnitude of their two-dimensional vector 
argument. Hence /^(a) = fsiWl) and Q{f3) = Q{\(3\). 

The formula ([1]) and the first term in formula ([2]) correspond to the naive picture 
(so-called optical approximation) implying that in the case of AA-coUision at the impact 
parameter (3 one can use the binomial distribution for N:^{f3) (see, for example, [211 122]): 

PoAN^) = C5 PiP f- , PiP) = 1 - Qif3) (T) 

with some averaged probability P{P) of inelastic interaction of a nucleon of the nucleus 
A with nucleons of the nucleus B. At that the P{f3) is considered to be the same for all 
nucleons of the nucleus A. In the optical approximation one has 

{N^mopt = AP{P) , V[N^iP)Ut = AP{P)QiP) . (8) 

The whole expression (jj]) for the variance is the result of more accurate calculation (see 
appendix R]) . when at first one calculates the probabilities of all binary NN-interactions, 
taking into account the impact parameter plane positions of nucleons in the nuclei A and 
B and only then averages over nucleon positions: 

V[N^m = {N^iPf) - {N^m' , (9) 

where 

„ B A 

(X) = {{X)b)a = XW TB{bk)dbk n TA{a,)da, . (10) 

fe=i j=i 

Here X is the average value of some variate X at fixed positions of all nucleons in the 
nuclei A and B; ( )^ and ( )b denote averaging over positions of these nucleons with 
corresponding nuclear profile functions. 

In the limit rjy <^ Ra, Rb the formulae ([5]) and Q reduce to 

fsia) ^ a^^TBia) , ^(^(ai, 02) ^ /(ai - as) ■ TB((ai + a2)/2) (11) 

with 

aj^j^ = Jdba{b) , I{a) = J dba{b) a{b + a) . (12) 

Note that in this limit the Q{f3) and hence the mean value ([T]) and the first term of the 
variance ([2]) depend only on the integral inelastic NN cross-section ajv^r, but the Q^^'^\(3) 
entering the second term of the variance ([2]) depends also on the shape of the function 
a{b) through the integral I (a) f|T2l) . 

Note also that using of the simple approximation with the 5-function: a (6) = (Xj^jySlb) 
for AW interaction gives the same result (as going to the limit rjy ^ Ra,Rb) only for 
the optical part of the answer, which is expressed through Q{P). If someone tries to use 
the approximation a{b) = (%Ar5(&) to calculate Q^^'^\(3), he will get I{a) = a^]yS{a) and 
9b = ct^nH'^i ~ '^2) ■ Tsiai), what leads to infinite Q^^'^\f3) at B > 2. Meanwhile, for 
any correct approximation of a{b) with a{b) < 1 (in correspondence with its probabilistic 
interpretation in classical Glauber approach) we find a finite answer for Q^^'^\/3). 
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Figure 1: The variance of the number of wounded nucleons in one nucleus for PbPb 
coUisions at SPS energies {aj^j^=31 mb) as a function of the impact parameter /3 (fm). The 
points • and ■- results of numerical calculations by the analytical formulae ([2])-(|l]), ( ITTj) 
and ( !T2|) using respectively the black disk (fT4|) and Gaussian ( ITSjl approximations for AW 
interaction; O and □ - results of independent MC simulations using for AW interaction the 
black disk ( IT^ or Gaussian (|T5l) approximation; * - the optical approximation result (E]) 
(the first term in formula ([2])); + - the Poisson variance: V[N:^{f3)] = {N^{f3)). The 
curves are shown to guide eyes. 

In the quantum case in Glauber approximation due to unitarity one has 

a{b) = (T^"(6) = a'"\b) -a'^{b) = 2lm-f{b) - \-f{b)\^ > , (13) 

where the 7(6) is the amplitude of AW elastic scattering. This leads to the restrictions: 
< cr*°*(6) < 4, < cr^'(6) < 4 and < cj^"(6) < 1. So in the quantum case the a{b) also 
admits a probabilistic interpretation |T9| [20]. 

In our numerical calculations we have used for a{b) the "black disk" approximation: 

a{b) = eir^ -\b\), (14) 

and Gauss approximation: 

a{b) = expi-byr-l) . (15) 

In both For the nuclear profile functions and Te we have used the 

standard Woods-Saxon approximation: 

TA{a) = Jdz p{r) , = + , p(r) = po (^1 + exp (16) 

with Ra = RqA^ , Ro=1.07 fm, k=0.545 fm and po fixed by the condition fdaTA{a) = l. 

The numerical evaluation of the contribution of the additional (contact) term in for- 
mula (E]) one can see in FigH] presented as an example for PbPb collisions at SPS energies 
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Figure 2: The same as in FigJH but for the mean number of wounded nucleons in one 
nucleus, calculated by formulae (II]), (E]), ([5]) and by independent MC simulations; * - the 
optical approximation result, calculated using formulae ([1]), (|3]) and (fT2l) . 

(rAr=lfm, o}YAr=31mb). For the control we have also carried out independent calculations 
of the mean values and the variances involved by MC simulations of the AA scattering 
presenting the results on the same figures. 

In Fig{T] we see that the contact term in ([2]) is essential and gives approximately the 
same contribution to the variance of the iV^(/3) in PbPb collisions at intermediate and 
large values of {3 as the first optical term. It's important that as we see in Fig{T]the results 
of independent MC simulations of the N^[(3) variance are in a good agreement with the 
results of the analytical calculations by formula ([2]) only if one takes into account the 
contact term. 

We see also in Fig{T]that for peripheral AA collisions at large /3, when P{P) becomes 
small (-P(/3) ^1, <5(/3) ~ 1), the optical approximation ([7]) reduces to the Poisson distri- 
bution with V[N^{(3)]^t ~ {N^m dH]). 

So only due to the contact term the variance of the N^[[3) is larger than the Poisson 
one for peripheral PbPb collisions (at /3 > 7 fm) in a correspondence with the indications, 
which one has from the experimental data on the dependence of multiplicity fluctuations 
on the centrality at SPS and RHIC energies [HH]. 

The week dependence of the results on the form of NN interaction at nucleon distances 
is also seen. In the case of using the black disk T\M approximation for a{h) the results lay 
systematically slightly higher, than in the case of using the Gaussian ( fT5l) approximation 
with the same value of a^^. 

In Figj2]we see that the mean value {N:^[[3)) ([1]), in contrast to the variance, coincides 
with the optical approximation result ([8]) and depends only on a^^ in the limit <^ 
RaiRb- The MC simulations also confirm this result. 

We would like to emphasize that the nontrivial term in the expression ([2]) for the 
variance arises only in the case of nucleus-nucleus collisions. At A = 1 or 5 = 1 it 
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Figure 3: An example of the loop diagram in AA-coUisions. 1 and 2 - nucleons of the 
nucleus A; 1' and 2' - nucleons of the nucleus B (see Ell 125] for details). 

vanishes. At A = 1 due to explicit factor A — 1 in ([2]) and at i? = 1 due to fact that 
in this case = Q'^{(3). This corresponds to the well known fact that for nucleus- 

nucleus collisions the Glauber approach doesn't reduce to the optical approximation even 
in the limit rjy ^ Ra, Rb (see, for example, [25]). 

The additional term, which arises in the expression for the variance in the case of 
nucleus-nucleus collisions, depends, as we have mentioned, not only on the integral value 
of inelastic NN cross-section = Jdba{b), but also on the shape of the function cr{b), 
i.e. on the details of NN interaction at nucleon distances, which are much smaller than 
the typical nuclear distances. In quantum Glauber approach it corresponds to the fact 
that in the case of AA collisions, in contrast with pA collisions, the loop diagrams of 
the type shown in Figj3] appear and one encounters the contact terms problem (see, for 
example, |23l[2l|25]). 

The second term in formula ([2]) is the manifestation of this problem at the classical 
level. In the case of a tree diagram the "lengths" of the interaction links in the transverse 
plane are independent. As a consequence the result expresses only through P(/3) - the 
probability of the interaction of a nucleon of the nucleus A with nucleons of the nucleus 
B averaged over its position in nucleus A. The P{(3) is the same for any nucleon of the 
nucleus A. In the case of the loop diagram in FigjSjthe "lengths" of the interaction links in 
the transverse plane are not independent and the result can't be expressed only through 
the averaged probability P(/3) and the correlation effects have to be taken into account. 

3 Variance of the total number of participants 

Now we pass to the calculation of the variance of the total number of participants 
V[N^{I3) + iy^(/3)] at a fixed value of the impact parameter f3. Clear, that for the 
mean value we have simply: 

and by iQ for the variance 

v[N^{P) + N^m = v[N^m + v[N^m + 2{{N^{msm - {N^miN^m} ■ 

(18) 

In naive optical approach there is no correlation between the numbers of participants 
in colliding nuclei at fixed value of the impact parameter: 
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Figure 4: The correlator between the numbers of wounded nucleons in colhding nuclei, 
calculated by analytical formulae (fT9|) - (!22|) and by independent MC simulations. The 
notations are the same as in FigJT] 



More accurate calculations fulfilled in accordance with (|9]) and (fTOj) (see appendix [B]) lead 
to 

{N^im^m - {N^miN^m = ab[q^''\(3) - , (i9) 

where 

Q(")(/3) = / dadbTA{a)TB{b)[l - fB{a+f3)f-'[l - fA{b-/3)f-'[l - a{a-b+/3)] , (20) 



Qi(3)= dbTB{b)[l-fA{b-P)\ 



and 



/a(6) = da TA{a)a{b- 



cr, 



NN 



TA{b) 



(21) 



f22) 



The Q{I3) and /^(a) are the same as in formulae ([3]), ([5]) and (ITT]) . Recall, that in our 
approximation /yi(6) = /yi(|^|) and Q{P) = Q{\f3\), then Q{/3) can be obtained from Q{/3) 
by a simple permutation of A and B. At A = B we have Q{f3) = Q{(3). 

The results of numerical calculations of the correlator 0191) by formulae fl2UI) - fl22p for 
PbPb collisions at SPS energies together with the results obtained by independent MC 
simulations of these collisions are presented in FigHl 

Comparing FigJH with Fig{T] we see that the contribution of the correlator to the 
variance of the total number of participants at intermediate values of /3 is about half of 
the variance for one nucleus V[N^{(3)] and is approximately equal to the contribution of 
the first optical term in ([2]). At large values of the impact parameter (/3 > 10 fm) the 
relative contribution of the correlator (fT9|) to the total variance (fT8l) is even greater. The 
results are again in a good agreement with the results obtained by MC simulations. (The 
small difference in the region 8-10 fm arises from the use of approximate formulae (fTTj) 
and 
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impact parameter, p (fm) 

Figure 5: The same as in FigJH but for the variance of the total number of wounded 
nucleons Nu,{f3) = N^{f3)+N^ {(3) in colhding nuclei. The variance V[Niu{l3)] is calculated 
by formulae ([2])-(j4]), (ITT]) . (fT2|) with taking into account the contribution of the correlator 
([T8])-([22D; + - the Poisson variance: V[N^{/3)] = {N^il3)). 

In FigslS] and [6] we present the final results for the variance of the total number of 
participants in PbPb collisions at SPS energies, taking into account the contribution of this 
correlator. (In Figj6]the same, as in Fig|5l but for the scaled variance: V[N^u{/3)]/ {Nyj{(3)), 
Nw{(3) = N^{f3) + N^{f3).) We see in particular that the calculated variance of the total 
number of participants is a few times larger than the Poisson one in the impact 

parameter region 8-12 fm. 



4 Variance of the number of binary collisions 

In this section we present the results of the calculation of the variance of the number of 
NN-collisions at a fixed value of the impact parameter /3 in the framework of the same 
classical Glauber approach [18] to nucleus-nucleus collisions. The details of calculations 
one can find in the appendix O 

As a result we found that the formula for the mean number of binary collisions again 
coincides with the well-known expression given by the optical approximation (compare 
with the formula fl29|) below): 

{N^oiim = ABx{(3) , (23) 

where 

X{(3) = J dadb TA{a)TB{h)a{a-h+P) ^ o-j^tv / TA{a)TB{a + P) (24) 

has the meaning of the averaged probability of NN-interaction. Numerically the mean 
value of the number of collisions as a function of the impact parameter /3 are shown in 

Figm 
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Figure 6: The same as in FiglHl but for the scaled variance V[N^{/3)]/ (A^^(/9)) of the total 
number of wounded nucleons in colliding nuclei, A''^(/3) = N^{f3) + N^{f3). 



In contrast to the mean value, the formula obtained for the variance of Ncoii{P)'- 
V[N,oiim = ABlxiP) + iB-l)xM + - (A+ fi- (25) 

differs from the optical approximation result (see below eq. (!30|) ). It depends not only on 
the x(/3) ([21]), but also on 



XiW) = J d<^ TA{a) J db TB{b)a{a-b+(]) ^ a^^ J da TA{a)Tl{a+P) 



(26) 



and 



Xi(/3) ^Jdb Tsib) 



da TA{a)a{a — b+/3) 



a^j^ J da TB{a)Tl{a+/3) . (27) 



The Xi is obtained from Xi by permutation of A and B. (Recall, that we consider the 
Ta and Tg depend only on the magnitude of their two-dimensional vector argument.) At 
A = B we have Xi = Xi- Note also that in the limit r^^t^RA, Rb the x, Xn X\ ^'^'^ hence 
the variance (l25ll depend only on cr^jy, but not on the form of the function (t(6) (it was 
not the case for the variance of the number of the wounded nucleons, see section [2] after 
the formula ffT2|l l 

For comparison we list below the optical approximation results, which assumes the 
binomial distribution for Ncoii{P) with the averaged probability x{P) of NN-interaction 
(see, for example, [2T| 122]): 



popt{N, 



coll I 



In this case one has 



x(/3)^-" [1 - xl/?)]^"""^-" 
{N^^mopt = ABx{(i) 



(28) 
(29) 
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Figure 7: The mean number of NN-collisions in PbPb interactions at SPS energies calcu- 
lated by the formulae ( 123|) and (124|) and by independent MC simulations as a function of 
the impact parameter /3 (fm). The notations are the same as in FigJT] 



and 

v[N,Mmopt = ABxm^ - xm = {NcouiPw - xm ■ m 

Note that for heavy nuclei xW) is small even for central collisions ~ ^nZ-^'a ^1); so 

the distribution fl28|) and the variance in optical approximation (l30l) practically coincide 
with the Poisson ones: V[Ncoii{l3)]opt ~ (NcoiiW)). 

Note also that in the case of pA interactions {A = 1 or B = 1) our result fl25|) for 
the variance of the number of collisions coincides with the formula ( l30l) obtained in the 
optical approximation. 

In Figs J8] and m as an illustration we present, the results of our numerical calculations 
of the variance of the number of collisions by analytical formulae (I24|) - (l27|) in the case of 
PbPb scattering at SPS energies together with the results obtained from our independent 
Monte-Carlo simulations of the scattering process. (In Figj9] the same as in FiglHl but for 
the scaled variance: V[Ncoui/3)]/ {Ncoii{f3)) .) 

We see that the calculated variance of the number of collisions at all values of the 
impact parameter /3 is a few times larger than the Poisson one, whereas the variance 
given by the optical approximation practically coincide with the Poisson one (see the 
remark after formula fl30|) ). The results obtained by independent Monte-Carlo simulations 
confirm our analytical result. (The small difference again can be explained by the use of 
approximate formulae (^^, fl2^ and 

We have also analyzed the dependence of the fluctuations on the diffuseness of the 
nucleon density distribution in nuclei. To study this dependence the calculations with a 
smaller (0.3 fm) than standard (0.545 fm) value of the Woods-Saxon parameter k (fT6|) 
were performed, what corresponds to the model of nucleus with a sharper edge (see FigsJTUl 
and [II]). 
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Figure 8: The variance of the number of NN-colhsions in PbPb interactions at SPS energies 
as a function of the impact parameter /3 (fm). The points • - results of calculations by 
analytical formulae (!24|) -( l27|) : * - the optical approximation result, calculated using 
formulae and (130|) : + - the Poisson variance: V[Ncoii{/3)] = {Ncoii{/3)). The notations 
are the same as in FiglH 



The calculations confirm that one would expect from simple physical considerations, 
more compact distribution of nucleons in nuclei does not change the mean number of 
wounded nucleons, but reduces its fluctuations, because in this case the number of 
wounded nucleons is more strictly determined by the collision geometry. As a result, 
the scaled variance of the number of wounded nucleons decrease with k, (compare the 
Figsl6]and[T0]). 

As for the number of binary NN-collisions, in this case due to more compact distribu- 
tion of nucleons in nuclei the mean number of collisions increases along with its variance. 
Therefore the scaled variance of the number of binary collisions weakly depends on the 
variation of the parameter k (compare the Figs 12] and Fig JTT]) . Important that in both 
cases the contribution of the contact term plays the crucial role. 

5 Discussion and conclusions 

It's shown that although the so-called optical approximation gives the correct results for 
the average number of wounded nucleons and binary collisions the corresponding variances 
can't be described within this approximation in the case of nucleus-nucleus interactions. 

In the framework of classical Glauber approach the analytical expression for the vari- 
ance of the number of participants (wounded nucleons) in AA collisions at a fixed value 
of the impact parameter is presented. It's shown, that along with the optical approxima- 
tion contribution depending only on the total inelastic NN cross-section, in the case of 
nucleus-nucleus collisions there is the additional contact term contribution, depending on 
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Figure 9: The same as in FigJHl but for the scaled variance V[Ncoii{P)]/ {Ncou{P)) of the 
number of NN-colhsions. 



details of NN interaction at nucleon distances. 

In classical Glauber approach this contact contribution arises due to taking into ac- 
count the interactions between two pairs of nucleons in colliding nuclei (a pair in one 
nucleus with a pair in another). It's found, that the interactions of higher order, than 
between two pairs of nucleons, don't contribute to the variance. Whereas the expression 
for the mean number of participants was proved to be exact already in the optical approx- 
imation, which bases on taking into account only the averaged probability of interaction 
between single nucleons in projectile and target nuclei. 

These results are obtained in the framework of pure classical (probabilistic) Glauber 
approach [18]. However it's possible to suppose, that in the quantum case the one-loop 
expression for the variance and the "tree" expression for the mean number of participants 
and binary collisions will be exact. 

Using obtained analytical formulae, the numerical calculation of the variance of the 
participants number in PbPb collisions at SPS energies was done as an example. Demon- 
strated that at intermediate and large impact parameter values the optical and contact 
term contributions are of the same order and their sum is in a good agreement with the 
results of independent MC simulations of this process. 

When calculating the variance of the total (in both nuclei) number of participants the 
correlation between the numbers of participants in colliding nuclei is taking into account. 
The analytical expression for the correlator at a fixed value of the impact parameter is 
obtained. The results of numerical calculations of the correlator for the same process 
of PbPb collisions show that at intermediate and large values of the impact parameter 
its contribution to the variance of the total number of participants is about half of the 
variance in one nucleus, again in good agreement with independent MC simulations. 

As a result for peripheral PbPb collisions the variance of the total number of partic- 
ipants, calculated with taking into account the contributions of this correlator and the 
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Figure 10: The scaled variance of the total number of wounded nucleons. The same as in 
Figjni but for the nucleon density distribution in nuclei ( IT6|) with a smaller value of the 
Woods-Saxon parameter k=0.3 fm. 



contact terms, occurs a few times larger than the Poisson one. 

In the framework of the same classical Glauber approach the analytical expression for 
the variance of the number of NN binary collisions in given centrality AA interactions 
is also found. Along with the optical approximation term it also contains other terms, 
which occur the dominant ones. 

Due to these additional terms the variance of the number of collisions at all values of 
the impact parameter is several times higher than the Poisson one, whereas the variance 
given by the optical approximation practically coincides with the Poisson one. Again the 
results obtained by the independent MC simulations confirm our analytical result. 

Important that these additional contact terms in the expressions for the variances 
arise only in the case of nucleus-nucleus collisions. In the case of proton-nucleus collisions 
they are missing and the variances are well described by the optical approximation. 

Note that we have used the simplest factorized approximation fl3Tl) for the nucleon 
density distribution in nuclei and do not take into account nucleon-nucleon correlations 
within one nucleus, which play a fundamental role, for example, in the description of 
particle production in nuclear collisions outside the domain kinematically available for a 
production from NN-scattering (so-called 'cumulative' phenomena) 

The additional contact contribution to the variance of the number of wounded nucle- 
ons, as we have found, arises due to interactions between two pairs of nucleons in colliding 
nuclei, which need to occur at the same position in the impact parameter plane. Taking 
into account nucleon-nucleon correlations within one nucleus must increase the proba- 
bility of such configurations and hence the contribution of the contact term. However, 
numerical accounting of these effects is beyond the scope of the present paper. 

Interestingly, the nontrivial contact terms in variances (missing in optical approxima- 
tion) arise in our approch already in the framework of the exploited factorized approxima- 
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tion for the nucleon density in nuclei, i. e. without taking into account nucleon-nucleon 
correlations within one nucleus. 

The authors thank M.A. Braun and G.A. Feofilov for useful discussions. The work 
was supported by the RFFI grant 09-02-01327-a. 



Appendices 

A Calculation of the variance of participants in one 
nucleus 

The geometry of y4i?-collision is depicted in FigJT2l All aj and are the two-dimensional 
vectors in the impact parameter plane. In the framework of the classical (probabilistic) 
approach [18] the dimensionless a{b) is the probability of inelastic interaction of two 
nucleons at the impact parameter value b (see also (fT3|) ). The Ta and Te are the profile 
functions of the colliding nuclei A and B. We are implying that for heavy nuclei the 
factorization takes place: 

A 

TAiai, ...,aA) = Y[TA{aj) . (31) 

i=i 

Convenient to introduce the abbreviated notation: 

J da = J TA^a) da = 1 . (32) 

All integrations imply the integration over two-dimensional vectors in the impact param- 
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Figure 12: Geometry of y4i?-collision. 



eter plane. In new notation the ( !T0|) takes the form 

(X) = ((X)^)^ = [xfldhU da, . (33) 

k=i j=i 

Recall that here X means average of some variate X at fixed positions of all nucleons in 
A and B; {)^ and ( )^ mean averaging over positions of these nucleons. 

We introduce the set of variates Xi, ...,Xa (each can be equal only to or 1) by the 
following way: Xj = 1, if j-th nucleon of the nucleus A interacts with some nucleons 
of the nucleus B and Xj = 0, if j-th nucleon doesn't interact with any nucleons of the 
nucleus B. The number of participants (wounded nucleons) in the nucleus A in a given 
collision at the impact parameter /3 is equal to the sum of these variates: 

N^W) = J:X,. (34) 
i=i 

Then we have for the mean value: 

A A 

(iV^(/5)) = E(^.) = E(TOB)A (35) 

and for the variance of N^{(3): 

V[N^m ^ (N^ipf) - {N^m' , {N^ipf) = ((E X,r) . (36) 

i=i 

At first we calculate the mean value (135|) . We denote by qj and pj the probabilities that 
the variate Xj will be equal to or 1 correspondingly. Clear that for given configurations 
of nucleons {aj} and {bk} in nuclei A and B: 

B 

^3 = 11(1- (Tjfe) , Pi = 1 - qj , (37) 

k=l 

where 

ajk = -bk + f3) (38) 
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and 

X- = 0-g, + l-p,=p, . (39) 
Note that pj and qj are the functions of aj, bi,...,bB and (3: 

q. = q.(aj, {6^}, /3) , pj = pj{aj, {fej, /3) . (40) 

Recall that we restrict our consideration by the region of the impact parameter /3 < 
Ra + Rb where the probability of inelastic nucleus-nucleus interaction cr4^(/3) is close to 
unity. Otherwise one has to introduce in formula ( 1371) for qj the factor l/cr4^(/3), where 

'^ABm = l-{{I[mi-cr,,))^)^ (41) 
j=ik=i 

and a^B = I ^(3 a^^i^P) is so-called production cross section, which can't be calculated in 
a closed form. 

Substituting (ETD-dSi) into ^ we have 

{N^m = A-j:{{q,)B)A- (42) 

Averaging at first on positions of the nucleons in the nucleus B, we find 

(g,)^ = (l-a,)^, 
where we have introduced the short notation: 

aj = I dbiaji = J dhiTB{hi)a{aj - h + (3) . (43) 

Averaging now on positions of the nucleons in the nucleus A, we have 

{{Qj)b)a= I da.il-a,)'' , (44) 



which is the same for any j, as the is the integration variable: 

{{^3)3) A = J da,TA{a,){l - a^f ^ Q(/3) . (45) 

Then by we find 

{N^m = A{l-Qm = AP{P) , (46) 
which coincides with formula ([1]) of the text, if one takes into account the connection 

^j = fBia,+P) (47) 

(see and (H5i) ). We see that the result for the mean number of participants (H^ is the 
same as in the optical approximation ([8]). 

We calculate now by the same way the variance of N^{f3). By fl36l) we have: 

{N^{Pf)= E {X,,X,,) + j:{X]) . (48) 
ii^i2=i i=i 
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Note that the (Xj-^Xj.^) can't be reduced to the product {Xj^){Xj^) . Just in this point 
the optical approximation breaks for AB coUisions. 
Since by dMD 

then for the first sum in fHHl) we find: 

A A 

E(^J) = E(^.) = ^N^m = mp) ■ (49) 

Because 

for the second sum in (HHi) using ( H5ll we have: 

Y: {X,,X,,) = AiA 2Q{P) + Q^''\P)] , (50) 

ii^i2=i 

where we have introduced 

Q'^""^^^)^ MA-l) ^ ^^^n<l32) b) A ■ (51) 

We calculate now Q'^^'^\(3). Averaging again at first on positions of the nucleons in 
the nucleus B, we have 

{q3.^32)B = {'^-^3.-^32+^^''''^)'' . 

where cxj^ and aj^ are given by P3l) and 

^{3132) ^ J dbiaj,iaj,i = J dbiTB{bi)a{aj, - 6i + (3)a{aj, - 6i + /3) . (52) 

Then averaging on positions of the nucleons in the nucleus A one can rewright (|5Tll as 
follows 

Q(i2)(/3) = J da,da2TAia,)TAia2)il - a, - + a^'^Y , (53) 

where by fl52|) 

= J rf\aiiff2i = y" c?6iTB(6i)or(ai-6i + /3)or(a2-&i + /3) = ^ij(ai+/3, a2 + /3) (54) 

(see notation ([6]) of the text). Substituting dH]), (gg]) and ([50]) into we find for the 
variance of N^((3): 

V[N^m = AQim - Qm + AiA-l)[Q^''\(3) - Q'm , 

which coincides with the formula ([2]) of the text if we take into account f j47l) . (l53ll and 
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B Correlation between the numbers of participants 
in colliding nuclei at fixed centrality 



The calculations are similar to ones in appendix |A] (we use the same notations). Along 
with the set of variates Xi, ...,Xa we introduce in the symmetric way the set of variates 
Xi,...,Xb (each can be again equal only to or 1). X^ = 0(1) if fc-th nucleon of the 
nucleus B doesn't interact (interacts) with nucleons of the nucleus A. Then similarly to 
for the number of participants (wounded nucleons) in a given event in the nucleus B 
we have: ^ 

N^W) = j:^k. (55) 

k=l 

Then 

(N^im^m = E E((^)b)a (56) 
j=i k=i 

and similarly to f l39|) 

XjXk = P,fc(l, 1) , (57) 

where the Pjk{^, 1) is the probability that the both variates Xj and X^ will be equal to 
1. For the probability Pjk{^, 1) one finds 

Pjkil, 1) = OTjk + (1 - crjk)pjkpjk , (58) 

where ajk is the probability of the interaction of the j-th nucleon of the nucleus A with 
the k-th nucleon of the nucleus B (see formula (138|) ) and pjk is the probability of the 
interaction of the j-th nucleon of the nucleus A with at least one nucleon of the nucleus 
B except the k-th nucleon (correspondingly pjj. is the probability of the interaction of the 
k-th nucleon of the nucleus B with at least one nucleon of the nucleus A except the j-th 
nucleon) : 

B A 
Pjk = 1-11 (1 - ^jk') , Pjk = 1 - Yl {1 - CTj'k) ■ (59) 

k>=i{k'^k) j'=Hj'^j) 

Combining ( I56l) -( l59l) and acting as in appendix |A] we find the formulae (fT9l) - (l22l) of 
the text. 



C Fluctuations of the number of collisions 

In this appendix we calculate the variance of the number of NN-coUisions in AB-interaction 
at fixed value of centrality in the framework of the approach under consideration. 

To calculate the number of collisions we define the set of the variates Yi, Ya, which 
can take on a value from to B. If in the given event the j-th nucleon of the nucleus A 
interacts with n nucleons of the nucleus B, then Yj = n. The number of NN-coUisions 
in the given event at the impact parameter f3 can be expressed through these variates as 
follows: ^ 

N,ouW) = T.Yj (60) 
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Clear that again (see appendix |A]) : 

P{Y, = 0) = q, = llil - a,k) (61) 

k=l 

To calculate P{Yj = n) for n = 1, i? we introduce {/ci, /c„} - the sampling from the 
set {1, ...,B} and {A;„+i, A;_b} - the rest after sampling. Then 

P{Yj = n)= Yl (^3k,--(rjk„{l - (rjk„+,)--{l - CTjk^) (62) 

{ki,...,kn} 

First we again calculate the mean value of the number of collisions: 

{Ncoum = j:mB)A- m 

i=i 

For a given configuration {aj} and {bk} we have: 

Y^ = j2nPiY, = n). (64) 

n=0 

Using (l62l) and averaging on positions of the nucleons in the nucleus B, one finds 

{yj)B = E ^ C'b^JII - c^.)''"" = Ba,. (65) 

n=0 

We use the same notations as in appendix El (see ( l43i) ). Averaging then on positions of 
the nucleons in the nucleus A, we finally find: 

{N,oum = ABxW) , (66) 

where 

x(/3) = J daioi = j daidbiffu = J daidbiTA{ai)TB{bi)a{ai — 6i + /3) , (67) 
and at r^r ^ Ra, Rb 

xiP) ^ (%Ar J daiTA{ai)TB{ai + (3) , (68) 



which coincides with the formulae ( 123|) and ( 12^ of the text. Comparing ( 166|) and ( l29l) 
we see that the result for the mean number of collisions is the same as in the optical 
approximation. 

In the rest of the appendix we calculate the variance of the number of collisions. To 
calculate the variance: 

V[NUP)] = {NLim - {NUP)f (69) 

one has to calculate 

{Nlum = iiJlY,)") = E {YnYn) + j:{Y') ■ (70) 
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So we have to calculate the following two sums: 

E {YnYn)= E {{YJJJ:)b)a (71) 

and 

E(>^') = E((>?)b)a- (72) 

To calculate the first sum we denote by k[, ...,k'^ - the indices of the nucleons of the 
nucleus B, which interact only with the nucleon ji of the nucleus A. By k'l,...,k'^ we 
denote the indices of the nucleons, which interact only with the nucleon of the nucleus A 
and by ki, ...,kr we denote the indices of the nucleons, which interact with both nucleons 
ji and j2- By ki, kB-n-m-r we denote the indices of the nucleons of the nucleus B, 
which don't interact with the nucleons ji and j2 of the nucleus A. Then the probability 
Pjj^j^ of such event in these notations is equal to 



where 



Pjlj2 — PjlPj2 5 (73) 
r n m B—r—m—n 

Ph = n ^iifc, n ^jiM 11(1 - aj.k^) n (^ - ^^nkj , (74) 



i=l i=l i=l i=l 

B—r—m—n 



Ph = n 11(1 - (Tj.fc/) n ^j2fc:' n (^ - ^i2fcj • (75) 

j=l «=1 j=l i=l 

Using (17^ and (175|l we can rewrite Pj^jg in the following form 

r n in B—r—m—n 

Phh = n ^iifc»^i2fc, n ^iife^(i-^i2fe^) n(i~^ji^f)^i2fcf n (i-f^jife,-o^i2fc.+f^jifc,c^j2fcj • 

i=l 1=1 «=1 1=1 

(76) 

The probability Pj^j^{n,m,r) that the nucleons ji and j2 of the nucleus A interact 
separately with n and m nucleons of the nucleus B and at that else simultaneously with 
r nucleons of the nucleus B is equal to 

Pjij-z in, m, r) = ^ pj^j^ , (77) 

where the sum means summing on all possible three sampling {k'^, ...,k'^}, {k'l, ...,k'^}, 
{fci, fcj.} from the set {1, ...,5}. After averaging (1771) on positions of the nucleons in 
the nucleus B we find 

(P,,,,(n,m,r))^ = , , ^' {y-zTix-zni-x-y+z)^-^--- , (78) 

n\m\r\[B — r — m — n)\ 

where we have used the short notations: 

X = a,-, , y = a,, , z = a^'^^'^ . (79) 

The (Tjj and aj^ are defined by ( H3i) and the a^^^^^^f is defined by ( l52ll in appendix El Then 
for the components of the first sum (I7T1) we have 

B B—r B—r—m 

(^OT^)b = EE E (m+r)(n+r)(P,,,,(n,m,r))5 . (80) 

r=0 m=0 n=0 
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After substitution of ( !78|) in ( IHOj) the lengthy but straightforward calculation leads to the 
simple answer 

{Y~Y~)s = Bz + B{B-l)xy = Ba^^'^^'> + B{B -l)aj^cTj^ . (81) 

For the components of the second sum (!72l) the similar but much more simple calculation 
gives 

{Y^)^ = Ba, + BiB-l)a] . (82) 
Averaging now on positions of the nucleons in the nucleus A, we can rewrite (170|) as 

{NLim = B 



B 



h¥'h=i i=i 
A{A-l) I daja2 ((a^^'^)^ + (5-l)(aia2)^) + A 1 rfai ((di)^ + (fi-l)(a^)^ 



Recalling now that (Xi, (J2 and cr'^^^^ are given by the formulae fH5]) and flM|) of the ap- 
pendix 1^1 we obtain 

{NUP)) = AB[x{/3) + {B-l)xM + {A-1)UP) + {A-l){B-l)x\f3)] (83) 

with x(/3), Xi(/3) and Xi(/3) defined by the formulae ([MD, (ES]) and ([27]) of the text. Using 
now the definition (169|1 and taking into account the formula (!66l) for {Ncoii{f3)) we come 
to the expression (l25l) of the text for the variance of the number of collisions. 
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